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ABSTRACT 

We have carried out spectroscopic observations in 4 cluster fields using Subaru's FO- 
CAS multi-slit spectrograph and obtained spectra for 103 bright disk field and cluster 
. galaxies at 0.06 ^ z ^ 1.20. Seventy-seven of these show emission lines, and 33 pro- 

OO ' vide reasonably-secure determinations of the galaxies' rotation velocity. The rotation 

velocities, luminosities, colours and emission-line properties of these galaxies are used 
to study the possible effects of the cluster environment on the star-formation history 
of the galaxies. Comparing the Tully-Fisher relations of cluster and field galaxies at 
similar reshifts we find no measurable difference in rest-frame i3-band luminosity at 
f~i i a given rotation velocity (the formal difference is 0.18 ± 0.33 mag). The colours of 

C^h' the cluster emission line galaxies are only marginally redder in rest-frame B — V (by 

0.06±0.04mag) than the field galaxies in our sample. Taken at face value, these results 
^ • seem to indicate that bright star-forming cluster spirals are similar to their field coun- 

I terparts in their star-formation properties. However, we find that the fraction of disk 

. galaxies with absorption- line spectra (i.e., with no current star formation) is larger in 

clusters than in the field by a factor of ^ 3-5. This suggests that the cluster envi- 
ronment has the overall effect of switching off star formation in (at least) some spiral 
. galaxies. To interpret these observational results, we carry out simulations of the pos- 

' sible effects of the cluster environment on the star-formation history of disk galaxies 

5^ , and thus their photometric and spectroscopic properties. This allows us to create mock 

samples of unperturbed "field" galaxies (with approximately constant star-formation 
rates) and perturbed "cluster" galaxies with different star-formation histories, includ- 
ing star formation truncation, with or without an associated starburst. We show that, 
if we select only bright galaxies with current star- formation (i.e,, with emission lines 
strong enough for rotation-curve measurements) , the average colours and luminosities 
of the "cluster" galaxies may not be very different from those of galaxies in the "field" 
sample, even though their star-formation histories may be significantly different. How- 
ever, the fraction of emission and absorption-line galaxies would change significantly. 
We also use these simulations to estimate the size of field and cluster galaxy samples 
that would allow us to differentiate the different star-formation scenarios considered. 
Finally, we find that the rest-frame absolute B-band magnitude of the field galaxies 
in our sample shows an evolution of —1.30 ± 1.04 mag per unit redshift at fixed ro- 
tation velocity. This indicates that the average SFR of bright disk galaxies evolves 
more slowly than the universal star-formation rate as determined from UV, Ha, far- 
infrared and radio studies. This suggests the evolution of the universal SFR density 
is not dominated by bright star-forming disk galaxies, in agreement with previous 
studies. 

Key words: galaxies: clusters: general - galaxies: evolution - galaxies: kinematics 
and dynamics - galaxies: spiral 
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1 INTRODUCTION 

Recent studies of distant galaxies have established that the 
fraction of SOs in rich galaxy clusters drops by a factor 2- 
3 from the local universe to z ^ 0.5, while the fraction of 
spiral galaxies increases comparatively (Couch et al. 1994; 
Dressier ct al. 1997; van Dokkum et al. 1998; Fasano et al. 
2000). Parallel studies have found a number of poststarburst 
galaxies within clusters (Dressier & Gunn 1983; Couch & 
Sharpies 1987), a suppressed star formation rate for spiral 
types relative to the field (Balogh et al. 1998), and evidence 
of an infalling spiral population out to z ~ 0.4 (Poggianti 
et al. 1999; Kodama & Bower 2001). This naturally leads to 
the scenario that cluster SOs could form from spiral galax- 
ies through interaction with the cluster environment (Jones, 
Sinail & Couch 2000; Kodama & Small 2001). 

Three mechanisms have been suggested for the trans- 
formation: galaxy-galaxy interactions, tidal forces, and 
gas stripping. Galaxy-galaxy interaction is most efficient 
in group environments, where relative velocities arc low 
(Zabludoir & Mulchaey 1998; Mulchaey & Zabludoff 1998; 
Ghigna et al. 1998), while tidal forces (Toomre & Toomre 
1972; Moore et al. 1996; Mihos, McGaugh, & de Blok 1997) 
are effective on small spirals. Gas stripping is expected to be 
an effective process in rich clusters (Abadi, Moore, & Bower 
1999; Quihs, Moore, & Bower 2000; Bekki, Couch, & Shioya 
2002; Vogt et al. 2004). Although each mechanism is effi- 
cient in different circumstances and/or environments, many 
of the rmmcrical studies of these mechanisms predict a star- 
burst could/should happen during the processes leading to 
the cessation of star formation required by the spiral-to-SO 
transformat ion . 

Until recently, however, there has been no compelling 
evidence indicating that spiral galaxies falling into dense 
environments experience a star-burst. Motivated by this 
Ziegler et al. (2003) investigated three clusters at 2 = 0.3- 
0.5 and used the TuUy-Fisher relation (TFR) to compare 
galaxy luminosities of field and cluster spirals at the same 
mciss, finding no measurable difference. In contrgist, Milvang- 
Jensen et al. (2003) carried out a similar study for the rich 
cluster MS1054.4-0321 at z = 0.83, and found evidence 
for some brightening in the cluster galaxies relative to the 
field ones at the same rotation velocity, implying more active 
star-formation of spiral galaxies in clusters. The situation is 
thus controversial, but these studies showed that using the 
TFR as a tool for comparing field and cluster galaxies at 
intermediate redshift could be successful. 

Following Milvang- Jensen et al., we have carried out a 
study of the TFR in nine additional distant rich clusters. 
Five have been observed with the VLT using the F0RS2 
multi-slit spectrograph (Seifcrt ct al. 2000) and the results 
(including a re-analysis of the MS1054.4— 0321 data) have 
been reported in Bamford et al. (2005A, B05A hereafter), 
where they have found further evidence for the bright- 
ening in the cluster spirals. In this paper, we report on 
our TFR study for the remaining four clusters observed 
with the Subaru Telescope using the FOCAS spectrograph 
(Kashikawa et al. 2002). Since our sample includes field 
galaxies, wc also discuss their evolution following Bamford, 
Aragon-Salamanca & Milvang- Jensen (2005B, B05B here- 
after). 

Before comparing the cluster spiral galaxies with the 
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field ones, one needs to be aware that the comparison could 
be complicated by strong selection effects. For example, the 
probability of observing a galaxy in a star-burst phase will 
be lower if the star-formation occurs over a short time-scale; 
on the other hand, a star-bursting galaxy will be brighter, 
and thus more likely to appear in a magnitude limited sam- 
ple. Similarly, one could think that if our sample contains 
galaxies in the fading phase after the starburst has ended, 
the average luminosity of our sample would be lowered, per- 
haps cancelling the effect of the initial brightening. However, 
it is unlikely that we will be able to measure rotation curves 
for galaxies after the cessation of star formation since emis- 
sion lines would be absent. All these effects will probably 
affect the average properties of the galaxies in our field and 
cluster samples, so we have carried out careful modeling for 
the correct interpretation of the observations. 

Throughout the paper, we assume flm ~ 0.3, Qa ~ 0.7, 
and Hq=7Q km s"^ Mpc'V AU magnitudes are in the Vega 
zero-point system. 



2 DATA 

2.1 Observations and data reduction 

Our sample contains rich clusters, mainly in the northern 
hemisphere, covering a wide redshift range, similar to B05A. 
Table lists the clusters observed, together with some of 
their properties and basic observational parameters. The 
MS2053.7-0449 cluster was also studied by B05A, so we 
should be able to test the consistency of the rotation veloc- 
ity measurements. 

Pre- imaging in the i?-band was carried out on 2002 June 
7 using FOCAS in imaging mode. The field-of-view was 6 
arcmin diameter, and the pixel scale was 0.1 arcsec pixel"^. 
Each field was taken at two different position angles to en- 
able flexible allocation of tilted slits (see below). The expo- 
sure times at each position angle were 180s for A2390 and 
240s for the other clusters. The seeing was 0.5-0.6 arcsec. 
The images at the different position angles were combined, 
and SExtractor version 2.0 (Berlin & Arnouts 1996) was 
used to identify the objects. The photometry was calibrated 
to Thuan & Gunn (1976) r-band (see SecIT^. 

For multi-object spectroscopy (MOS) targeting, we de- 
termined the priority of each galaxy according to how many 
of the following conditions were met: (i) disky morphology 
in pre-imaging (must), (ii) bright enough to observe, (iii) in- 
clination > 45° , (iv) HST images available from the archive, 
(v) magnitudes and/or the redshift available from literature, 
and (vi) known to be a cluster member. We designed the 
masks by allocating slits to galaxies with the highest pri- 
orities. The slits were to be placed along the major axis 
of galaxies to measure rotation velocities. Two masks with 
complementary position angles were designed for each clus- 
ter except MS2053.7— 0449 for which we decided to use only 
one mask. This allowed us to use slit angles < 45° relative to 
the spatial direction for most galaxies, avoiding degrading 
the resolution beyond acceptable limits. A few slits included 
more than one galaxy by chance. Table Ogives the details 
of the masks. In total 108 slits were placed on 116 galaxies 
in 7 masks, of which 4 galaxies were observed in two masks 
and 112 were unique. In the following we always refer to the 



number of unique galaxies unless stated. The morphological 
distribution of the selected galaxies is discussed in Sec. 12.51 
The procedures followed for target selection, priority allo- 
cation and mask design were very similar to those used by 
B05A. 

The spectroscopic observations were carried out on 2002 
Aug 10-11 using FOCAS in MOS mode. The 300B grism was 
used, yielding a dispersion of 1.4Apixel~^. The slit width 
was 0.6 arcsec along the dispersion axis for all the objects, 
achieving a spectral resolution R ~ 1200, slightly higher 
than B05A. The lower limit of the spectral range was set 
to 4700 A by the order blocking filter, and the upper limit 
was up to 9400 A, depending on the geometrical position of 
each slit on the mask. The seeing, determined from the spec- 
tra of ~ 2 stars observed in each mask, was 0.4-0.7 arcsec, 
and we binned x2 along the spatial direction, achieving a 
final spatial sampling of 0.2 arcsec pixel"^ for untilted slits. 
Given the differences in redshift, we made shorter exposures 
for A2390 and longer ones for MS2053. 7-0449 to achieve 
roughly the same depth in every cluster. Table Ogives the 
exposure time of each mask. 

Bias subtraction, flat flelding, connection of the two 
CCD detectors, and distortion correction were done us- 
ing the FOCASRED package, developed by the instru- 
ment team. Wavelength calibration and sky-subtraction 
were carried out with standard IRAF tasks to obtain the 
2-dimensional spectra. 



2.2 Redshift and cluster membership 

The centre of each galaxy on the 2D spectrum was identified 
by the peak of the continuum component without emission 
lines. Central spectra were extracted for each galaxies in 0.2- 
0.6 arcsec apertures, and the redshift was measured using 
emission and/or absorption features. When more than one 
emission line were well detected, the mean of the redshifts 
was taken. The emission lines included [O II]A3727, 115, H7, 
H/3, [O III]A4959, [O III]A5007, [N II]A6548, Hq, [N II]A6583, 
[S II]A6716, and [S II]A6731. The redshifts of 77 galaxies 
were measured in this way. When emission features were not 
seen, absorption features were used instead. We employed a 
cross-correlation method with a model template spectrum. 
Redshifts for 27 absorption-line galaxies were measured in 
this way. We could not determine the redshifts of 8 galaxies 
because of poor detection of their absorption features or 
dubious emission line identification. We drop these galaxies 
from further analysis, and the remaining is 104. The redshift 
range of the galaxies was 0.06 ^ 2 ^ 1.20 with median of 
z = 0.42. When considering only the galaxies in the final 
TFR sample (see Sec. 12. 4L the redshift range was 0.06 5^ 
z ^ 0.74, with a median of 2 = 0.39. 

The field-of-view of FOCAS was too small to cover the 
whole of the clusters, so cluster membership was decided on 
the basis of velocity only, without considering the spatial 
location of the galaxies. Galaxies with velocities within ±3cr 
from the velocity centre of the cluster were classified as clus- 
ter members. The fraction of cluster and field galaxies were 
48% and 52%, respectively. The numbers for each cluster are 
summarised in Table |5| 
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Table 1. Summary of the clusters and spectroscopic observations. 



Cluster 




R.A. 


Dec. 


z 


cr 


Mask 


P.A. 


T 

^ cxp 


Seeing 












(J2000) 


(J2000) 




(km s~^) 


ID 


(deg) 


(s) 


(arcsec) 








A2390" 




21 53 36.8 


-1-17 41 32 


0.2282 


1294 


PI 


0° 


1200 X 4 


0.7 


12 


14 


1 














P2 


-45° 


1200 X 3 


0.5 


13 


13 




MS1621.5H 


-2640" 


16 23 34.5 


-1-26 34 17 


0.4271 


735 


PI 


45° 


1800 X 4 


0.4 


15 


17 


1 














P2 


-45° 


1800 X 4 


0.6 


15 


16 




MS0015.9H 


-1609" 


00 18 33.5 


-1-16 26 03 


0.5490 


984 


PI 


0° 


1800 X 5 


0.5 


16 


16 


2 














P2 


90° 


1800 X 5 


0.5 


20 


22 




MS2053.7- 


-0449'' 


20 53 44.6 


-04 49 16 


0.583 


817 


PI 


0° 


1800 X 5 


0.4 


17 


18 





" Position, z and cr from Girardi & Mczzctti (2001). 

^ Position and z from Stocke et al. (1991), a from Hoekstra et al. (2002). 
° Number of galaxies observed in two masks. 



2.3 Absolute B-band magnitudes and inclinations 

To supplement the pro-imaging photometry and obtain 
colour information, wo obtained magnitudes for the targeted 
galaxies from the literature. We also collected WFPC2 im- 
ages from the archive, as in B05A. The available sources 
of photometry, other than our own FOCAS prc-imaging, 
were the following: for A2390, the CNOC survey (Yec ct al. 
1996) and WFPC2 images {F555W, F702W, FSUW); for 
MS1621.5-I-2640, CNOC (EUingson et al. 1997) plus WFPC2 
{F555W, F814W^); for MS0015.9+1609, CNOC (Ellingson 
et al. 1998) plus data from Dressier & Gunn (1992), and 
WFPC2 (F555W, F8141^); for MS2053.7-0449, WFPC2 
images (F606VF, F702W, FSUW). Note that not aU these 
data were available for all of our targeted galaxies (see be- 
low). 

The CNOC survey used the Thuan & Gunn photomet- 
ric system and provides r-band total magnitudes and g — r 
colours within 6.4 arcsec apertures. Dressier & Gunn (1992) 
used Schneider, Gunn, & Hoossol (1983) photometric sys- 
tem, and give total (/-band magnitudes and g — r and r — i 
colours within the fitting radius. The WFPC2 images were 
extracted from the axchive with the basic reductions com- 
pleted. We ran SExtractor for these images and obtained 
AUTOMAGs, which provide good estimates of total mag- 
nitudes within 0.1 mag. This uncertainty in the total mag- 
nitude is too small to have any effect on our TFR analysis. 
The zero-points were calibrated with reference to the lat- 
est WFPC2 manual. The available HST photometric bands 
differed from cluster to cluster and object to object, but 
at least two bands were always available, so we had some 
colour information for every galaxy within the WFPC2 field. 
The redder bands were used to define the aperture centres 
when deriving colours. The FOCAS pre-iniaging was not 
photometrically-calibrated at the telescope, so that photom- 
etry w£is calibrated using the gala^xies that were in common 
with the CNOC. A constant offset was used to transform 
SExtractor's AUTOMAG in the prc-imaging to the total 
Thuan & Gunn r-band magnitudes from CNOC. The error 
of this calibration was ~ 0.08 mag, estimated from the r.m.s. 
of the transformation. Thus, the magnitudes derived from 
the FOCAS pre-imaging are in the Thuan & Gunn r-band 
system. 

For each galaxy, the magnitude in the band that was 
closest to the rest-frame Johnson-Morgan B-band was cho- 
sen among those available as the starting point to calculate 



rest-frame B magnitudes. To obtain a representative colour 
we used the data from the CNOC, WFPC2, and Dressier 
& Gunn (1992) (in that order of priority). The Galactic ex- 
tinction of Schlegel, Finkbeiner, & Davis (1998) was applied 
to these magnitudes and colours. The magnitude was then 
converted to the absolute rest-frame B-band magnitude by 
referiug to the galaxy colours and K-corrections of Fukugita, 
Shimasaku, & Ichikawa (1995). The concordance cosmology 
(cf. §1) was used to calculate the distance modulus. In many 
cases we had several colours available from different sources, 
allowing us to estimate that the error in the rest-frame B- 
band magnitude associated with using different colours was 
0.03 mag on average. Moreover, we estimate that the error 
associated with using different photometric bands to esti- 
mate rest-frame _B-band was 0.08 mag on average. This in- 
dicates that the values of our rest-frame B magnitudes are 
quite robust, with uncertainties of the order of 0.1 mag or 
less. 

In a few cases the only available photometry came from 
the FOCAS pre-imaging. In that case, the rest-frame B mag- 
nitude was estimated from the pre-imaging magnitude ap- 
plying an Sc galaxy K-correction. Nine objects belong to this 
category, with only four malting it to the final TFR study. 

The inclination of the disk component was measured 
using GIM2D (Simard et al. 2002) following B05A. The 
WFPC2 images were used whenever available, while the prc- 
imaging was used otherwise. When the WFPC2 images were 
used, we always selected the band with the highest S/N, 
yielding the highest confidence. There were 48 galaxies with 
an inclination i > 40° (i.e., large enough for reliable ro- 
tation velocity determination) which had their inclinations 
measured in both the pre-imaging and the WFPC2 F814W 
images. For these galaxies, the average ratio of the (sini)~^ 
values determined from the pre-imaging and HST data re- 
spectively was 1.02 ± 0.01. Similarly, there were 39 galaxies 
with i > 40° that had their inclinations determined from 
F555W and F814W. For these, the average (sini)~^ ratio 
was 0.98 ± 0.01. Hence the systematic error in mixing the 
inclinations obtained from the different WFPC2 bands and 
the pre-imaging is only ^ 2%, and can be safely ignored. 
Note that the FOCAS images had very good seeing, and the 
value of the inclination is dominated by the outer isophotes, 
making our ground-based data perfectly adequate for incli- 
nation determination at these moderate redshifts. 

Finally, the internal reddening was corrected as a func- 
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Figure 1. Absolute B-band magnitude vs. redshift for the 104 
galaxies with measured redshifts. The different symbols indicate 
absorption-line galaxies (crosses), emission-line galaxies (trian- 
gles), and emission-line galaxies with measured rotation velocity 
(circles). Small and large symbols denote field and cluster galax- 
ies, respectively. The two marks with arrows are the galaxies at 
z = 1.195 (cross) and 1.041 (triangle). The dashed line a,t z = 0.19 
is used to define an unbiased sample for our analysis (see text for 
details) . 

tion of inclination following the prescription of TuUy & 
Fouque (1985). The absolute extinction at face-on position 
was assumed to be 0.27 mag in B following B05A. The de- 
rived absolute B-band magnitudes (Mb) are shown in Fig.Q 
as a function of redshift. 

There were seven galaxies in common with the VLT 
sample of B05A in the MS2053. 7-0449 field. Our B-band 
magnitudes for these agreed well with those by B05A, with 
five out of the seven being within 0.1 mag. This ensures the 
consistency of the derived magnitudes even though we used 
partly different procedures and photometric sets from B05A 
in our derivation. 

Our sample does not contain enough bright galaxies at 
the lower end of the redshift range due the small observed 
volume. Since this could cause a bias in our analysis, we 
will only consider galaxies with z > 0.19 in our study. In 
doing so we ensure that our field and cluster samples contain 
galaxies covering similar luminosity ranges (—22.5 < Ms < 
— 19.5). The redshift and magnitude ranges studied here are 
comparable to the 'matched' sample of B05A (z > 0.25 and 
Mb < —19.5), although we do not reach the redshift of their 
most distant cluster (MS1054.4-0321 at z = 0.83). 

2.4 Rotation velocities 

In order to be able to measure the rotation velocities of the 
77 galaxies with detected emission lines (see Table I^J , we 
first fitted the continuum component and removed it in each 
spectrum by interpolation. Then every emission line was cut 



out to a postage stamp image, and input to the software 
ELFIT2PY (B05A) together with the inclination and the ab- 
solute B-band magnitude to measure the rotation velocity. 
ELFIT2PY is based on the algorithm of ELFIT2D (Simard 
& Pritchet 1998, 1999), and calculates the best fitting the 
rotation velocity and scale-length of the emission assuming 
an intrinsic rotation curve and an exponential hght distri- 
bution. The software takes into account the instrumental 
resolution and the seeing, as well as the slit width (opti- 
cal beam smearing). The output rotation velocity has to be 
corrected for inclination. The inclination is also necessary 
as an input to calculate the light distribution through the 
slit. The absolute B-band magnitude is taken into account 
in the fitting process since we assumed the universal rota- 
tion curve of Persic & Salucci (1991). The output rotation 
velocity would change, on average, by ~ 10 km s~^ if we 
adopted a flat rotation curve, the exact value depending on 
the luminosity /mass of the galaxy (cf. B05A). Such small 
change would not affect the conclusions of this study. The 
spatial centre was determined from the centroid of the con- 
tinuum, while we allowed ±1 pixel flexibility in the velocity 
centre (~ 76 km s"^ at 5500 A and z = 0) to let ELFIT2PY 
account for possible errors in redshift and wavelength cali- 
bration. In fact, of the objects used in the flnal TFR (see 
below) 85% had offsets of 0.2 pixels or less. 

As part of the rotation velocity determination we car- 
ried out a quality-control process to ensure that only reliable 
measurements were used in the analysis. Our process is sim- 
ilar to the one used by B05A, so we will only describe here 
briefly the main quality-control steps. First, we did not at- 
tempt to measure the rotation velocities when the emission 
was nuclear (i.e., spatially-unresolved) or when the mean 
S/N was below a given threshold. Second, when measure- 
ment yielded a rotation velocity smaller than its estimated 
1(7 error, we rejected the measurement as unreliable. Third, 
the fits were visually inspected and those clearly wrong were 
also rejected. The main causes for rejection included: not 
reaching the fiat rotation or the turn-over point (cf. Ver- 
heijen 2001), clear non-exponential distribution of the emis- 
sion, strong asymmetry, and a deep stellar absorption close 
to the emission. In a few cases the slits appeared to have 
been cut in the wrong direction, so these galaxies were also 
rejected. Note that our treatment of the turn-over point is 
different from that of B05A, and we will discuss this issue 
in Sec. 13] After this process, only 33 galaxies (43%) with de- 
tected emission lines yielded secure rotation velocities that 
we can reliably use in the final TFR analysis. For reference, 
we attempted to fit 63 galaxies using ELFIT2PY, of which 
15 were rejected due to clear velocity under-estimates by 
ELFIT2PY or no sign of turn-over, 9 for low quality fits, 2 
for lack of rotation structure, and 4 for problems not directly 
related to the rotation velocity measurement. 

When more than one emission line were securely mea- 
sured for a galaxy, we used the velocity measured using 
the line with the highest S/N, while B05A took the error- 
weighted mean value. The number of emission lines mea- 
sured per galaxy was 2.5 on average for the 33 galaxies, and 
25 of these had more than one emission line. The ratio of 
the rotation velocities determined with the two highest S/N 
lines has a mean of 0.97 and a standard deviation of 0.11. 
Hence the rotation velocities of the different emission lines 
were consistent with each other typically within 11%. 
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In Fig.|5|we show the observed rotation curves of the 33 
galaxies in the final sample, together with the fitted models, 
in order to give a visual impression of the quality of the data 
and the fits. We note that a few galaxies may have irregular 
kinamatics (e.g. C0016_P1 11_A). Examples of objects re- 
jected by our quality-control procedure are shown in Fig. |21 
The rejected objects were not used in our analysis. 

The measured rotation velocities were corrected for in- 
clination multiplying them by (sini)"^. The uncertainties 
in the final values were estimated taking into account the 
errors in the inclinations from GIM2D, and the errors in the 
rotation velocities from ELFIT2PY. Table |21 provides some 
statistics for the spectroscopic data. 

The internal consistency of the rotation velocities was 
checked using four galaxies observed in two different masks. 
Three of them showed emission lines. After the quality- 
control process, two galaxies were left with only one accept- 
able spectrum due to slit angles being more than 60 degrees 
and/or having rotation curves of insufficient quality from the 
other spectrum. Thus, for only one galaxy we measured two 
good-quality rotation velocities from two independent spec- 
tra (M1621_P1 12_A and M1621_P2 10_A). These rotation 
velocities were 163tl^ and 186^^4 (kms~^) respectively (be- 
fore inclination correction). These measurements were thus 
consistent (~ 1.3a discrepancy), and the second value was 
chosen for the final sample due to marginally better S/N in 
the spectum. 

The excellent seeing of the Subaru observations and the 
relatively high spectral resolution of the FOCAS data al- 
lowed us to test the effect of the seeing and the instrumen- 
tal resolution in the measurement of the rotation velocities. 
This is important when comparing with the VLT results (cf. 
B05A), obtained in poorer seeing and with slightly lower 
spectral resolution. The spectra of the surviving 33 galax- 
ies were convolved with Gaussians to simulate a ~ 1.0 arc- 
sec seeing and a spectral resolution R ~ 950. These values 
roughly correspond to the VLT data of B05A. The rotation 
velocities were then measured in the same way. The ratio of 
the rotation velocity determined from the convolved spec- 
trum to the original one was 1.02 ± 0.03, and no systematic 
trend was found as function of fitted emission scale- length. 
Hence we should be able to compare our data with that of 
B05A without any correction. We note that this applies only 
to the galaxies with secure rotation velocities in our sample, 
and does not necessarily mean that the method we have used 
to measure rotation velocities is not affected by the seeing 
conditions. 

Among the seven galaxies that we have in common with 
the VLT data of B05A fSec. 12.31 . three galaxies yielded se- 
cure rotation velocities. Unfortunately, only one of these, 
M2053_P1 07_A in Fig.H was measured securely by B05A. 
Our velocity estimate for this galaxy was 182t;? km s'^ 
before inclination correction, while B05A obtained 154^jJ 
km s~^. Although this is our only direct comparison, the re- 
sults are consistent with each other within 1.4a. The fully- 
corrected magnitudes, rotation velocities, and colours of our 
sample are summarized in Table |2l 



2.5 Morphologies 

The morphology of the galaxies was double-checked on the 
HST images after the observations were completed. Al- 
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Figure 3. Examples of rotation velocity measurements which 
have been identified as insecure and thus rejected. The reasons 
are written at the left-upper part of each panel. The symbols and 
lines are the same as in Fig. |^ 



Table 4. Visual Morphology of 62 galaxies for which HST images 
are available. 





Cluster 


Field 










Total 


Total 


15 13 


34 


62 


SO 


<1 <6 


<2 


<9 



though the FOCAS imaging allowed us to reliably select 
disk galaxies, it didn't have enough resolution to separate 
SOs from spirals. Among the 104 galaxies targeted, 62 had 
been imaged with HST. Some statistics of our visual inspec- 
tion are given in Table 0] As we targeted disk galaxies (see 
Sec. 12. IL our sample includes no elliptical galaxies except 
for one that happened to be on the slit of another target and 
was identified as a cluster member of A2390. Sometimes it 
was difficult to distinguish SOs from spirals even with the 
HST-images. Although we did not find any clear SO galax- 
ies, we did find some galaxies for which the classification 
was dubious. For that reason, Table 3] gives an upper limit 
to the number of SOs in our sample. By design, most galaxies 
in our sample turned out to be spirals. We estimate that at 
most ~50% of the galaxies with absorption-line spectra were 
SOs. We will come back to this issue later. We note that all 
the galaxies in the final TFR sample (i.e., with secure rota- 
tion velocities) that had HST images had clear spiral mor- 
phologies, as expected. Postage stamps of the HST/FOCAS 
images of the galaxies in the final sample are presented in 
Fig. El 
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Figure 2. The observed rotation curves of the 33 galaxies that survive our quality-control process and are used in the final TFR 
analysis. The cluster and field galaxies are displayed in the first two and the last three columns, respectively, in redshift order. The 
velocities are "observed", as opposed to "intrinsic", i.e., not corrected for inclination, non-zero slit width, seeing effects, etc. The 
vertical dash line shows the spatial centre. The large circles denote the velocity centre determined from Gaussian fits to each spatial 
column with the peak higher than 2 times the background RMS and the FWHM between 0.7 and 1.4 the resolution. The crosses 
are the same as the circles but with the peak as low as 0.7 times the background RMS and the FWHM between 0.5 and 2.0. Dotted 
lines show the ELFIT2PY fitted model when the peak of the Gaussian is above 0.01 times the background RMS. The zero-point of 
the velocity in each panel is set to the adopted redshift for the object. The small circles and horizontal ticks at the bottom of each 
panel show the flux distribution of the spectrum and model fit, respectively. The Flux is derived by integrating under the fitted 
Gaussians. One pixel corresponds to 0.2(cos6)~^ arcsec where 6 is the slit tilt angle. The width of the box is 12.2 arcsec. Labels 
at the top-left of each panel indicate whether the galaxies belong to the clusters or the field. The numbers right below these labels 
indicate the residuals, in Mg, from the local TFR (cf. Fig.|S] see Sec. l3.lt . The cluster name, mask ID, slit ID, redshift, and psr 
fitting element of the galaxies are indicated on the middle-left of each panel. The emission line ID is shown on the top-right corner. 
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Table 2. Summary of the emission line measurements. 



Observed Emission- ^ rot- 

detected measured 



Cluster 












N/ 
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A2390 




16 
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MS1621.5H 


-2640 


11 


18 
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5 


17 


3 


10 


MS0015.9H 


-1609 


17 


16 


3 


11 


15 


5 


5 


MS2053.7- 


-0449 


6 


12 





4 


11 


1 


4 


Sub-total 
Total 




50 


54 
112 


8 


28 


49 

77 


13 


20 

33 



Number of galaxies in clusters. 
^ Number of galaxies in field. 

Number of galaxies with unknown z. 




log [km/s] 

Figure 5. The TFR of the 33 galaxies for which reliable rotation 
velocities were measured. Filled and open circles are the cluster 
and field galaxies at 2 > 0.19. Field galaxies at 2 < 0.19 (crosses) 
are kept separated because of possible sampling bias (see Sec. 12. 31 . 
The dotted line shows the local TFR of PT92. The solid and 
the dashed lines show the linear fits to the cluster and the field 
galaxies at z > 0.19 with the slope fixed to that of PT92. 

3 RESULTS 

The primary purpose of this study is to see whether spiral 
galaxies in clusters have higher star formation activity than 
those in the field. In the following, we focus on three different 
aspects to investigate it: magnitude as function of rotation 
velocity (i.e., TFR), colour, and fraction of spiral galaxies 
with absorption- line spectra (i.e., no current star formation). 

3.1 Tully-Fisher relation 

The TFR of our sample is shown in Fig. |3 There are 6 
field galaxies ai z < 0.19. Although the number is small. 



they are consistent with the local TFR of Pierce & TuUy 
(1992; PT92 hereafter; Mb = -7.48 log Kot - 3.10 - 0.27^). 
This suggests that our low redshift data is consistent with 
the results of PT92. The galaxies at z > 0.19 tend to be 
brighter than the local TFR of PT92 when compared at the 
same rotation velocity. This suggests a possible luminosity 
evolution at higher redshifts. More on this later. 

No clear difference is apparent between the cluster and 
the field galaxies at z > 0.19. Using a simple mini- 
mization method with internal errors taken into account (cf. 
B05B), we fitted linear regression lines with the same slope 
as the local TFR of PT92 to the field and cluster galaxies. 
We find that the cluster and the field galaxies are 0.64 ±0.24 
and 0.46 ± 0.23 mag brighter, respectively, than the local 
spirals in PT92. The cluster galaxies appear to be 0.18 mag 
brighter than the field galaxies, but the difference has no 
statistical significance. 

An obvious problem with the above analysis is that our 
TFR has been derived for galaxies at different redshifts. If 
there is luminosity evolution with redshift, it should be taken 
into account. To assess this issue, in Fig.|S|we plot the resid- 
uals (in magnitude) of our TFR from the local one of PT92 
as function of redshift. There is a trend suggesting positive 
luminosity evolution with redshift. A linear fit to the field 
galaxies at z > 0.19 gives 

AMs = (-1.30 ±1.04) 2 -I- 0.09 ±0.46. (1) 

The median redshift of the sample used is 2; = 0.39. The 
fitted line has AMb ~ at z = 0, consistent with the lo- 
cal TFR of PT92. The six field galaxies at z < 0.19 are 
also consistent with the fit. Even taking into account the 
possible luminosity evolution with redshift. Fig. |S] also sug- 
gests that there is no obvious difference in the behaviour 
of the cluster galaxies from that of the field galaxies. The 
intrinsic standard deviations around the TFR linear fit, af- 
ter the luminosity evolution with redshift is taken off, are 
0.73 mag and 0.78 mag for the cluster and the field galax- 

1 The intercept of the TFR in PT92 may be 0.45 ± 0.12 mag 
brighter if the difference between Hq that was derived from the 
sample of PT92 (Pierce 1994, 86 km s-^ Mpc"!) and the one 
used in this paper is taken into account (B05B). Our conclusions 
are not affected by this issue. 
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Table 3. The data for the 104 galaxies in our sample with measured redshift. The galaxies are sorted by redshift first for the cluster 
members and then for the field galaxies. The columns are: (1) cluster, mask and slit IDs; (2) R.A. and (3) Dec; (4) redshift; (5) 
cluster (C) or field (F) membership fiag; (6) sample status ('TFR'=in the final TFR; 'em'= emission-line galaxy without secure 
rotation velocity; 'abs'= absorption-line galaxy); (7) inclination (90° =edge-on); (8) absolute rest-frame _B-band magnitude, (9) 
rotation velocity, and (10) rest-frame B — V colour. 
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0.56 

0.88 
0.40 
0.93 
0.79 
0.96 
0.84 

0.90 
0.65 

0.77 
0.96 
0.50 
0.96 
0.59 
0.40 
0.82 
0.55 
0.71 
0.69 
0.82 
0.70 
0.40 
0.90 
0.55 
0.50 
0.73 
0.71 
0.84 
0.64 
0.86 
0.67 



ies, respectively. Hence, there is also little difference in the 
scatter of the TFR of cluster and field galaxies. For our sam- 
ple, this scatter is about twice that of the local galaxies in 
PT92. This is probably due to the fact that local TFR stud- 



ies concentrate on samples of nice, regular spirals, while the 
high-redshift samples are, perforce, less clean. 

In Fig. Q we show the cumulative distribution of our 
TFR residuals from the local relation of PT92. We exclude 
the galaxies at z < 0.19 to avoid possible biases (although 
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Table 3. Continued. 



Cluster/Mask/Slit R.A. Dec. z Mem. TFR i Mg logVrot B-V 

[J2000] [J2000] [deg] [mag] [km s"!] [mag] 
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PI 
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0.71 
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0.67 


C0016 


P2 


11.A 


00 


18 


20.1 


+16 


26 


00 


0.1670 


F 


TFR 


77 


-17.86 + 0.07 
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0.46 


M1621 


PI 


10^ 


16 


23 


34.8 


+26 


35 


45 


0.1721 


F 


em. 


71 


-17.97 + 0.05 
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0.2459 


F 


em. 


63 


-21.68 + 0.05 
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0.4471 


F 


em. 


72 


—21.26 ± 0.07 
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the results would not change if we included them). The dis- 
tribution of the cluster galaxies resembles that of the field 
galaxies. The same conclusion is observed after the measured 
luminosity evolution with redshift is taken off by Eq. Q . A 
Kolmogorov-Smirnov test indicates that the significance of 



the difference between the cluster and the field galaxies is 
56% for the raw values, reducing to 21% after the evolution 
is subtracted. 

One might suspect that in the process of rotation ve- 
locity determination and quality-control (Sec. 12.41 we may 
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Figure 4. HST/FOCAS images of the galaxies in our final TFR sample. We show HST images (labeled with band IDs) when 
available, and FOCAS images (marked 'preim') otherwise. The size of the box is 12". 2 X 12". 2, the same as the spatial dimension 
used in Fig.|3 The panels are arranged as in Fig. 1^ 



have missed a number of cluster galaxies that are undergo- 
ing a strong star-formation episode with irregular kinemat- 
ics. A simple test for this possible effect is to compare the 
fraction of emission-line galaxies which have secure rotation 
velocities for the cluster and the field samples. According to 



Table|21 the respective fractions for cluster and field galaxies 
are 46 ± 13% and 41 ± 9%, using Poissonian errors (see also 
Ziegler et al. 2003). Since these fractions are very similar, 
there is no clear indication that such a selection bias is in- 
troduced by our quality control. In conclusion, there is no 
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0.4 0.6 
Redshift 

Figure 6. The TFR residuals from the local relation of PT92 as 
function of redshift. The symbols are the same as in Fig.|S] The 
dashed line shows a linear fit to the field galaxies at z > 0.19. 
The solid line is the same but for the cluster galaxies with the 
slope fixed to that of the field galaxies. 




-2 2 

AMg [mag] 



Figure 7. The cumulative distribution of the TFR residuals from 
the local relation of PT92 (2 > 0.19 galaxies only). The solid 
and the dashed lines in panel (a) show the cluster and the field 
galaxies, respectively. The number at the left-top indicates the 
significance of the difference between the cluster and the field 
distributions from a Kolmogorov-Smirnov test. Panel (b) is the 
same as panel (a) but the measured luminosity evolution with 
redshift is subtracted using Eq. 0. The right panel also shows 
the cluster (circles) and field (crosses) galaxies from the 'matched' 
sample of B05A after subtracting the evolution in the same way 
(see text for details). 



statistically-significant difference in the TFR between the 
cluster and field galaxies of our sample. 



3.2 Rest- frame B — V colours 

We have seen that the TFR of the cluster galaxies in our 
sample looks similar to that of the field galaxies. In the previ- 
ous section we found no evidence for a selection bias against 
star-burst galaxies introduced by our quality-control pro- 
cess. A further test for possible selection biases can be car- 
ried out if we look at the colour distributions of our samples, 
since galaxies with a strong star-formation episode would be 
very blue. Fig.|H|shows the rest-frame B—V colour (Johnson- 
Morgan system) for all the galaxies in our sample with 
redshift and colour information available. The rest-frame 
colours were derived as described in Sec. The resuhs 
below do not change if we use observed colours which are 
independent of K-correction effects. In the figure, the clus- 
ter line-emission galaxies with rotation velocities unavailable 
(large triangles) do not show bluer colour compared with 
those with rotation velocities available (large circles). The 
mean rest-frame B — V colour of these galaxies are 0.63±0.05 
and 0.61 ± 0.04, respectively. Thus, the galaxies dropped in 
the rotation velocity measurement process are not signifi- 
cantly bluer than the ones we kept. Hence it is unlikely that 
we have statistically dropped cluster galaxies with higher 
star- formation from the rotation velocity measurement. 

The possible existence of such bias for the field galax- 
ies used in our TFR can be tested in the same way. The 
mean rest-frame B — V colour of the field galaxies in the 
TFR is 0.55 ±0.03 mag, while that of the field emission-line 
galaxies without secure rotation velocities is 0.57±0.03 mag. 



Thus the field galaxies used in the TFR sample well the field 
emission-line galaxies in terms of their colour distribution. 

The colours allow us a further consistency check of our 
TFR result obtained in Sec. 13. II Higher star-formation rates 
would produce bluer colours at the same time as excess B- 
band luminosity. The average rest-frame B — V colour of 
the cluster galaxies used in the TFR is 0.61 ± 0.04, while for 
the field galaxies it is — 0.55 ± 0.03. These values are only 
1.2a apart. If the cluster galaxies had higher star-formation 
than the field ones (as suggested by B05A), we would ex- 
pect them to be bluer, and not redder. Furthermore, Fig. El 
shows no correlation between AMb and (B — V)rcst for ei- 
ther the field or the cluster galaxies. Such a correlation could 
be expected if excess star-formation were responsible for any 
excess luminosity in these galaxies. The distribution of the 
cluster and field galaxies in this diagram is indistinguish- 
able. These two points are consistent with the fact that the 
TFR of the cluster galaxies in our sample is similar to that 
of the field galaxies. 

As a final check, if we include all the emission galaxies 
in our colour test (i.e., not only those that make it into the 
TFR analysis) we still do not find any significant difference 
between the cluster and the field. The colour of all the emis- 
sion galaxies in the clusters (0.62 ± 0.03) is slightly redder 
than that of field galaxies (0.56 ± 0.02), but the difference 
is again not significant. We will address quantitatively the 
colour and magnitude changes expected in different evolu- 
tionary scenarios in Sec. 15.11 

3.3 Spiral galaxies with absorption-line spectra 

The presence of emission lines is a direct diagnostic for cur- 
rent star formation. The fraction of field spiral galaxies with 
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0.4 
Redshift 

Figure 8. Rest-frame B — V colour distribution of our sample 
galaxies. We plotted 93 of the 104 galaxies with known redshift, 
after dropping 9 galaxies with no colour information (see Sec. lLj.ijt 
and two galaxies at z > 0.8. The symbols are the same as Fig. Ill 
The dotted lines show the colours of E, SO, Sab, Sbc, Scd, and 
Im types at 2 = from Fukugita et al. 



no emission lines detected is 9 ±4% (using Poissonian errors; 
see Tablel^J. For all the cluster galaxies (including suspected 
SOs), the fraction is much higher, 44±9%. Even if we exclude 
all the suspected SO galaxies and the serendipitous elliptical 
(50% of the cluster galaxies with no emissions. Sec. 12.51 , we 
estimate that in our sample at least 29 ± 9% of the spiral 
galaxies in the clusters have no emission lines. Thus, there is 
an excess of non star-forming spiral galaxies in the clusters 
compared with the field. Our data suggest a factor of at least 
~ 3 difference, but it could be as high as a factor of ~ 5. 
The existence of such 'passive' spiral galaxies in rich clus- 
ters is consistent with the results of Dressier et al. (1999): 
if (at least some) spirals transform into SOs when they fall 
into clusters, their star formation (and therefore their emis- 
sion lines) must be switched ofi' at some point. This result, 
coupled with the lack of difference found in the TFRs of the 
cluster and field spirals in our sample (Sec. 13. ll . seems to 
suggest that the process leading to the eventual cessation 
of star-formation is not accompanied by a drastic star-burst 
event. This point will be explored further in Sec. 15.11 

The majority of the cluster absorption-line galaxies 
have colours typical of early-type galaxies, and the peak 
of their colour distribution is clearly different from that of 
emission-line galaxies. The mean rest-frame B ~ V colour 
of the absorption-line galaxies in our clusters is 0.76 ± 0.04, 
which is 0.20 ± 0.05 mag redder than the average colour of 
the field emission-line galaxies. A model star-burst contain- 
ing 20% of the total stellar mass and solar metallicity reaches 
{B — V)rest ~ 0.7 about 1 Gyr after the burst ended, and 
evolves very slowly thereafter ((-B — V^)rcst ~ 0.8 at ~ 2 Gyr; 
Bruzual & Chariot 2003). A truncated star- formation model 
behaves similarly. Thus the absorption-line galaxies in our 
sample have typically not experienced any significant star- 
formation for ~ 1 Gyr. However, there are a few absorption- 
line galaxies in our clusters that show bluer colours and are 
thus expected to be younger. This indicates that the sam- 
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Figure 9. The TFR residuals vs. rest-frame B — V colours for the 
33 galaxies with secure rotation velocities. The redshift evolution 
has been subtracted using Eq. 0. The zero-point of the colours 
is set to the mean of the field galaxies. The symbols are the same 
as Fig. 1^ The solid and dotted lines show the results of the sim- 
ulation discussed in Sec. 15. if indicating the regions where 68% of 
the cluster emission-line galaxies are distributed in models 9 and 
1 respectively (see text for details). The dashed line is the same 
but for the field galaxies. 



pie of absoption-line disk galaxies in our clusters contains 
galaxies at different stages in their evolution. 

If we limit the above analysis to those galaxies with 
bona-fide spiral morphologies as determined from the HST 
images, the mean B—V colour of the absorption- line galaxies 
in our clusters becomes 0.58 ±0.10. This is closer to the field 
emission-line galaxies, and would suggest a time sequence for 
the morphological transformation from spirals to SOs: spirals 
first switch off their star formation (i.e., no emission lines), 
and then become redder as their stellar populations age and 
their morphologies are transformed. 



4 COMPARISON WITH PREVIOUS WORK 
4.1 Detailed comparison 

Before we can consider the implications of our results it is 
important to assess how reliable the rotation velocity mea- 
surements of intermediate redshift galaxies are. To carry out 
an external check, we compare our results with those of other 
work which uses a similar technique at comparable redshifts. 
We first compare the TFR of our field galaxies with the ones 
derived by other workers in Fig. 1101 We plot the TFR resid- 
ual vs. redshift to separate redshift evolution from other 
effects. We do not compare our results with studies which 
use emission line widths (e.g., Rix et al. 1997) instead of re- 
solved rotation curves because that method is substantially 
different from ours, potentially causing significant system- 
atic differences. 
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Vogt et al. (1996, 1997) investigated tiie TFR of 16 field 
galaxies at 0.15 < 2 < 1 selected based on 7-band magni- 
tudes from Forbes et al. (1996), the Deep Extragalactic Evo- 
lutionary Probe project (Mould 1993), and the Hubble Deep 
Field (Williams et al. 1996). Eleven of the rotation curves 
obtained were labeled 'high quality' by them. The observed 
position-velocity diagrams were published for all the galaxies 
of their sample, and our visual inspection criteria also class 
these 11 as good-quality rotation curves. A visible turn-over 
point in the rotation curve was one of their quality-control 
criteria, as it is for us. Their results are consistent with the 
ones we found for our field galaxies. 

Simard & Pritchet (1998) studied 12 field galaxies with 
0.26 ^ 2 ^ 0.43. The sample was selected by its strong [O II] 
emission (equivalent widths > 20 A). The B-band absolute 
magnitude range of their sample was similar to that of Vogt 
et al. (1997). The measured rotation velocities were not cor- 
rected for inclination, so we applied an average correction to 
their galaxies before placing them on Fig. 1101 by assuming 
(Vobs) = 0.7854(Vodgc-on), which holds when galaxies are 
randomly oriented (cf. Simard & Pritchet 1998). We also re- 
move two galaxies for which the derived error is larger than 
the derived rotation velocity. Since these authors did not 
publish the observed position-velocity diagrams, we could 
not compare our visual quality-control with theirs. We do 
not know whether a visible turn-over point of the rotation 
curves was one of their quality-control conditions. The TFR 
residuals of their sample are, on average, slightly brighter 
than in our field sample, although their redshift range and 
their sample size are small. 

Ziegler et al. (2003) investigated the TFR of 13 cluster 
and 7 field galaxies with 0.3 ^ z ^ 0.6 in three cluster fields. 
All the observed position-velocity diagrams were shown in 
Jager et al. (2004), and 16 out of their 20 rotation curves 
pass our quality-control criteria, while four are rejected as 
rotation velocity overestimates. These authors also required 
a visible turn-over point in the rotation curves as one of their 
acceptance criteria. The TFR residuals of their field galaxies 
are consistent with the local galaxies of PT92, even at inter- 
mediate redshifts. Therefore, their field galaxies are slightly 
fainter than ours at a given rotation velocity, although their 
sample size is small. 

Bohm et al. (2004) investigated the TFR of 77 field 
galaxies with redshifts between 0.1 and 1 in the FORS Deep 
Field. This work includes the galaxies studied by Ziegler et 
al. (2002). The sample was selected in the R band, requir- 
ing that the spectral energy distribution of the galaxies was 
later than E/SO. They selected 36 of their 77 galaxies as high 
quality, and for 18 of them the observed position- velocity di- 
agram was shown. We agree with their quality assessment 
for these 18 galaxies, which are probably their best. A visible 
turn-over point of the rotation curve was one of their con- 
ditions for high quality, while those galaxies with a smaller 
extent or asymmetries were classified as low quality. The 
TFR residuals of their sample is consistent with that of our 
field sample. Naturally, our Eq. Q agrees well with their 
resulting fit [see their Eq. (11)]. 

B05B, using the same sample as B05A (which includes 
also the data published by Milvang- Jensen et al. 2003), stud- 
ied the properties of field spirals at intermediate redshift. 
Their sample contains 89 field galaxies with 0.1 ^ z < 1. 
The observed position- velocity diagrams of six galaxies were 



shown as an example in B05A, and our quality-control ac- 
cepts five of them, while one is rejected as an underestimate. 
The main difference in the quality-control process of B05A 
with ours is that they placed less weight on the condition 
of a visible turn-over point, because ELFIT2PY should take 
that into account when estimating the errors (see B05A for 
details). For reference, we rejected 57% of the galaxies in 
our sample during the quality control process, while B05A 
rejected 47% of theirs. The TFR residuals from B05B are, 
at all redshifts, about 0.5 mag brighter than those from our 
sample, while the rate of redshift evolution (i.e., the slope 
in Eq. is consistent. On panel (b) of Fig. |7|we compare 
the distribution of the TFR residuals from our Subaru sam- 
ple with the 'matched' sample of B05A after removing the 
redshift evolution in the same way. A Kolmogorov-Smirnov 
test indicates that the distributions for field galaxies from 
both samples disagree with 98% confidence. 

To finalise our comparison with previous work, we con- 
sider TFR studies that compare cluster and field galaxies. 
Ziegler et al. (2003) found no difference in the TFR between 
cluster and field galaxies, consistent with our result. We note 
that the behaviour of both the field and cluster galaxies of 
Ziegler et al. (2003) in the TFR residuals vs. redshift dia- 
gram is rather similar to that of ours. On the other hand, 
B05A (including Milvang- Jensen 2003 data) , found that the 
average rest-frame _B-band luminosity of cluster spirals was 
significantly (~ 0.7 mag) brighter than that of field ones at 
the same rotation velocity. We do not find that effect in our 
Subaru data. Panel (b) of Fig. |7| shows a comparison of our 
sample with that of B05A. 

4.2 More on quality control issues 

The visual quality control process is, perforce, subjective. 
The main difference found with some of the previous stud- 
ies discussed above is the treatment of the turn-over point of 
the rotation curve. We find reasonable consistency between 
our results and those from studies that took the visible pres- 
ence of a turn-over point in the observed velocity-position 
diagram as a requirement. It is therefore reasonable to ex- 
pect that differences between our results and those from 
other studies (e.g., B05A) could be due to this issue. 

Requiring a visible turn-over point has both advantages 
and disadvantages from a scientific point of view. If we re- 
quire it, we may drop from our sample a number of small 
galaxies which undergo strong star formation (cf. 'downsiz- 
ing' scenario) , or cluster galaxies which may have more con- 
centrated star-formation due to interactions with the cluster 
environment (cf. Moss & Whittle 2000). On the other hand, 
if we estimate the rotation velocity of a galaxy using a region 
that is too small (without reaching the turn-over point), the 
dark matter may not yet dominate, thus underestimating 
it. (e.g., Verheijen 2001). ELFIT2PY is designed to over- 
come this problem by direct 2-dimensional model spectrum 
fitting. So far, however, the possible systematics introduced 
by the software when no visible turn-over point is visible 
have not been explored. This is the main reason why, in this 
paper, we decided to use the visibility of the turn-over point 
as one of our quality-control conditions. 

Detailed model simulations to test the behaviour of 
ELFIT2PY under different conditions are being carried out, 
and will be discussed elsewhere (Bamford et al. in prepara- 
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tion). Nevertheless, we have found that galaxies with small 
fitted emission-line scale-length tend to have brighter TFR 
residuals. Fig. Illf a) shows the Subaru TFR residuals as a 
function of the ratio between the fitted emission-line scale- 
length and the photometric disk scale-length. We show both 
the galaxies with secure rotation velocity determinations 
and those with insecure ones (i.e., the ones that did not 
pass our quality-control) . The photometric disk scale-length 
is derived using GIM2D on the redest HST images, when 
available, and pre-imaging otherwise. If we look at all the 
emission-line galaxies in the figure, there is a trend in the 
sense that objects with smaller emission scale-lengths have 
brighter offsets in the TFR residuals. A possible explana- 
tion for this trend is that ELFIT2PY may tend to un- 
derestimate the rotation velocity for galaxies with smaller 
emission scale-lengths, but it could well be that these differ- 
ences are real. On the other hand, the galaxies with secure 
rotation velocities (i.e., the ones used in the analysis) do 
not show such a trend, so this gives us some confidence in 
our quality-control/rejection process. Nevertheless, if con- 
centrated emission is linked with increased star star forma- 
tion (and thus brightening), eliminating galaxies with com- 
pact emission from our analysis could hide some interesting 
and important physical effects. It is clear that using the pres- 
ence or absence of a turn-over point in the rotation curve 
as a quality-control criterion is not as clear cut as one could 
naively expect, and it could significantly affect the conclu- 
sions. 

In our sample, a significant fraction of the galaxies with 
a fitted emission scale-length smaller than the photomet- 
ric disk scale-length tend to have insecure rotation velocity 
measurements. Many of them were rejected because they 
clearly underestimated the rotation velocity and/or have 
strong nuclear emission (cf. H2.41 . It could be argued that an 
acceptance threshold should be introduced in the emission- 
line scale-length. A reasonable threshold could be set by 
0.8 X (photometric disk scale-length) if our quality-control 
can be taken as optimal (cf. Fig. lllb ). However, we decided 
not to impose such a strict limit and to use visual inspection 
to reject/accept measurements because we may lose some 
galaxies whose rotation velocities can be determined reason- 
ably well but which have more centrally-concentrated star 
formation. These galaxies could provide important clues of 
the evolutionary processes taking place in the cluster envi- 
ronment (Moss & Whittle 2000; Milvang- Jensen 2003). This 
issue will be discussed in a future paper. 

Fig. Illf b'l shows an apparent correlation between the 
TFR residuals and the rotation velocity. We note that 
this correlation would be the natural consequence of the 
TFR residuals being a function of the rotation velocity 
[AMb = 7.48 log Vrot + ...]. The arrow on the panel shows 
the direction of the correlated errors. Although Bohm et 
al. claimed the same trend at log K ot > 1-8 as evidence of 
mass-dependent evolution (using a bootstrap bisector fit), 
the trend can also be explained by random errors in the 
rotation velocities for a magnitude-limited sample without 
requiring any evolution (see B05B). The insecure rotation 
velocity measurements tend to occupy the lower velocity 
region of this trend. This may be an indication that the 
rotation velocities of these galaxies could be systematically 
underestimated, although is is not possible to rule out that 
a real effect is present. 
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Figure 10. The TFR residuals of our field galaxies vs. redshift 
compred with those obtained from other studies. The cosmol- 
ogy and the face-on _B-band extinction corrections are adjusted 
to match the ones adopted in this paper. The different sym- 
bols indicate: field galaxies in our sample (filled circles); high- 
quality galaxies from Vogt et al. (1996, 1997) (triangles); Simard 
& Pritchet (1998) (open stars); field galaxies from Ziegler et al. 
(2003) (asterisks); high-quality galaxies of Bohm et al. (2004) 
(open circles); and field galaxies from B05B (crosses). The dashed 
line is the same as in Fig. |S] 



To summarise the above discussion, the subjectivity 
of the quality-control process, especially the treatment of 
the turn-over point of the rotation curve, can significantly 
change the TFR results, and it could be the source of dis- 
crepancies between different studies. It is thus important to 
present the data, including the observed position-velocity 
diagrams, in such a way that cross-checks between different 
samples and studies can be carried out. Of course, we must 
not forget other possible sources of discrepant results such as 
different instrument/telescope combinations or real physical 
differences between samples. 

The discussion here illustrates the difficulties in com- 
paring datasets observed by different groups and/or using 
different instruments. Even with the utmost care, inconsis- 
tencies are often unavoidable. But even if we cannot compare 
the galaxies in our Subaru sample with the ones in the VLT 
sample of B05A, the differences between field and cluster 
galaxies within each dataset should be more reliable because 
their field and cluster samples were selected, observed and 
analysed in exactly the same way. Since we cannot rule out 
that the differences between cluster and field galaxies found 
by B05A may be real, at this stage we will take their re- 
sults at face value, as a working hypothesis, and explore the 
implications. 
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Figure 11. Panel (a): the TFR residuals as a function of the 
the ratio between the fitted emission-line scale-length (i?cm) and 
the photometric disk scale-length (i?disk) for a-H th^ galaxies with 
detected emission lines in the Subaru sample. The redshift evo- 
lution was subtracted using Eq. Q. Thick large and thin small 
circles denote the cluster and field galaxies with secure rotation 
velocities, respectively. The solid and open triangles represent the 
cluster and field galaxies with unreliable (i.e., rejected) measure- 
ments. The dotted line indicates where Rem = 0.8 x -Rdisk- Be- 
low this line most measurements have been rejected (triangles). 
Panel (b): the TFR residuals as function of rotation velocity. The 
correction for inclination in the velocity is taken off to see the 
raw measured values. The arrow shows how an error in the ro- 
tation velocity changes the position of galaxies in this diagram. 
The direction and the dimension of the arrow corresponds to a 
36% underestimate in the rotation velocity at a fixed magnitude, 
which results in making the TFR residuals 1 mag brighter. 



5 IMPLICATIONS OF OUR RESULTS 

We have seen no obvious difference between the cluster and 
the field emission galaxies in the TFR. This result is consis- 
tent with the results of Ziegler et al. (2003). We have also 
found no difference in the rest-frame B — V colour distri- 
butions of the field and cluster emission-line spirals. On the 
other hand, we have detected a difference in the fraction of 
absorption- line spirals in our field and cluster samples. This 
is rather a puzzling situation, since this last result is usu- 
ally interpreted as evidence for a transformation of spiral 
galaxies into SOs in cluster environments, while the first two 
results seem to suggest no difference in the cluster and field 
emission-line spirals. It is not impossible that, for a spiral 
galaxy sample like the one we have studied, the statisti- 



cal effect on the colours and luminosities produced by the 
galaxy transformations we have postulated is too weak to 
be detected. This issue can only be addressed with detailed 
modeling. 

The situation seems even more complicated if we accept 
the results of the study published by B05A, who found that 
cluster spirals are brighter than the field ones at the same 
rotation velocity, suggestive of a period of enhanced star- 
formation before it is switched off. Why don't we detect 
the same phenomenon? One possibility is that cluster-to- 
cluster differences could play an important role, although 
both studies concentrate on relatively rich clusters. Another 
possibility is that, due to the small number statistics, our 
Subaru sample may have missed several cluster galaxies that 
are undergoing this putative starburst, preventing us from 
detecting the same effect found by BQ5A. The likelihood of 
such possibility needs to be investigated. 

To explore these two issues, we have computed several 
simulations of the expected changes in the statistical proper- 
ties of field and cluster galaxy samples under different plau- 
sible evolutionary scenarios. 

5.1 Model expectations 

The basic assumption of our models is that field spiral galax- 
ies fall into clusters and get transformed into SOs after their 
star formation is extinguished by the cluster environment 
(e.g., after their gas is removed by ram-pressure stripping). 
We start by assuming that the parent field galaxy is forming 
stars at a relatively constant rate. The star-formation his- 
tory of such galaxy after it enters the cluster will depend on 
the details of its interaction with the intra-cluster medium, 
the cluster tidal field, other galaxies, etc. One possibility is a 
complete and sudden truncation of the star formation, as ex- 
pected if all the gas is removed from the disk and halo of the 
galaxy. A second possibility is that the star formation ceases 
gradually (e.g., exponential decay), as expected if only the 
halo gas is removed while the gas disk is retained. In this 
case, the star formation will decay as the available gas is 
used up and not replenished from the halo reservoir. A third 
possibility is that in the interaction process, the halo gas is 
removed, and the disk gas is compressed, producing a period 
of enhanced star-formation leading to a rapid exhaustion of 
the gas reservoir. Our models explore these possibilities, and 
predict the observational consequences we would expect in 
a study like ours. 

For the field galaxies we assume stellar populations with 
constant star formation and a fixed age of 5 Gyr. The re- 
sults are not sensitive to this value provided that is larger 
than a few Gyrs. For the cluster galaxies we use the same 
model until the galaxy encounters the cluster. At that time, 
different star formation histories are considered, with differ- 
ent star-burst contributions and timescales. Fig. 1121 shows 
some examples of the star- formation histories we used. The 
star formation after the encounter with the cluster has an 
e- folding timescale r, with burst mass fraction fb relative 
to the underlying population. Our calculations assume solar 
metallicity for the stellar populations, and use Bruzual & 
Chariot's (2003) models with a Chabrier (2003) initial mass 
function and the Padova 1994 stellar evolutionary tracks. 
Table |^ lists the values of the adopted model parameters 
(see below). 
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Figure 12. Examples of star-formation histories explored in 
our simulations. Model IDs correspond to those in Table ISl The 
absolute scale of the vertical axis is arbitrary. 



We further assume that the parent field galaxies have a 
rest-frame B-band luminosity function like that determined 
from the 2dF survey by Norberg et al. (2002) for the local 
universe (a = -1.21 and Mg = -20.43 with h = 0.7), and 
make Mg evolve with redshift according to Eq.Q, while 
keeping a constant. We select galaxies randomly from this 
luminosity function to generate a field population and a clus- 
ter population whose star-formation histories come from the 
models described above. We assume magnitude uncertain- 
ties of 0.1 mag, similar to the ones in our data. The mag- 
nitude limit of the model samples is set by the observed R- 
band magnitudes to mimic our observational selection. At 
each redshift, this i?-band limit is translated into a value of 
Mb (z) using our adopted cosmology and assuming the spec- 
tral energy distribution of the field galaxy models. Because 
of the way our observed samples were selected, this limit is 
reasonably close to Alg for the 'quiescent' field population. 
Obviously, we take into account that the star-formation his- 
tory of a galaxy alters its luminosity and colours, bringing 
in or out of our simulated sample depending on its apparent 
magnitude. 

We also assume that the field spirals retain their spiral 
morphology and emission line spectrum during the whole pe- 
riod covered by our simulations (i.e., they continue forming 
stars at a constant rate). This assumption seems reasonable 
for our sample since 91% of our field galaxies have detected 
emission lines (cf. Sec. 13.311 . In the case of the cluster galax- 
ies, we assume that the spiral morphology is observable for 
a time Atsp after their entrance in the cluster/start of the 
possible burst. After that time, the cluster galaxies do not 
enter our sample since they would not be classified as spi- 
rals anymore. That defines the maximum age of the cluster 
galaxies in our simulated samples. It is important to remem- 
ber that the absorption- line galaxies in our cluster sample 
could be contaminated by SOs when considering the results 
of our simulations. 

Thus, a population of infalling cluster spirals is built 
in our simulations by randomly-sampling the field luminos- 
ity function and assigning to the galaxies an age tt counted 
from their entrance in the cluster/start of the burst. This 



age is taken randomly from the interval (0 ^ tt, ^ Atsp). 
The corresponding magnitude change in each band is deter- 
mined from the model. Obviously, if tt is comparable to or 
smaller than r, the galaxy would be brighter than the par- 
ent galaxy because of the starburst. At later times, it would 
become fainter due to the cessation of star formation and 
the aging of the stellar population. The magnitude limit de- 
termines whether the galaxy enters the sample or not. If the 
galaxy ends up in the sample, a rotation velocity is assigned 
to it using the TFR of PT92 assuming luminosity evolution 
from Eq.Q. In this process the luminosity that we use to 
determine the galaxy's rotation velocity is the one it had 
before it entered the cluster. The underlying assumption is 
that the interaction with the cluster environment alters the 
galaxies' luminosities (via a change of their star formation 
histories), but not their masses/rotation velocities. The er- 
ror in logVrot (km s~^) is set to 0.1, comparable with our 
data. A comparison sample of field spiral galaxies is built 
by randomly populating the field luminosity function down 
to our apparent magnitude limits, and assigning rotation 
velocities to them in a similar way. 

We further define the parameter Atern (0 ^ Atem ^ 
Atsp) by requiring that emission lines are observables only 
if tb ^ Atem- This parameter controls the fraction of 
absorption- line spiral galaxies. Note that this is not a free 
parameter since the equivalent width (EW) of the emission 
lines (and therefore their detectability) depends on the rel- 
ative intensity of the star formation at any given time, and 
therefore it is possible to estimate reasonable values of Atem 
for each model (see below). 

The observed redshift Zobs is set to 0.4, which is close 
to the median of our sample (Sec. 12.21 . We set Atsp = 3 
Gyr, based on the numerical simulations of gas stripping by 
Bekki, Couch, & Shioya (2002), and r = 0.1 Gyr following 
Quilis et al. (2000). For each model, Atem is determined 
by considering reasonable observational limits in the EW 
of the emission lines. The number of ionising (Lyman con- 
tinuum) photons, NLy, is taken from the population syn- 
thesis models, and used to calculate the Ha luminosity as 
L{Ha) = 1.36 X 10"'^'^ • NLy (erg s~^). Our calculations as- 
sume that all the ionising photons from the hot stars are 
absorbed (case B recombination), an electron temperature 
Te = 10'' K and an electron density Ue = 100 cm~'^. We 
further apply a mean extinction of Ay = 1.1 mag to the 
emission lines, as derived for nearby SDSS spirals by Naka- 
mura et al. (2004). Seaton (1979) extinction law is then used 
to derive the EW of the line from the line and continuum lu- 
minosities of the models. The observational threshold is set 
to lOA for Hq, which is more or less similar to that of our 
observations. Note that because the extremely rapid evo- 
lution of the EW after the star-formation has ceased (see, 
e.g., Alonso-Herrero et al. 1996), a factor two difference in 
the EW limit does not affect Atem significantly. Equally, 
setting the limit using a different line would make very little 
different. For the burst strength we use a reasonably high 
value of 0.2 initially. The model calculated with this set of 
parameters is called our 'standard' model (number 1). We 
will explore later the effect of varying these parameters (cf. 
Table I^J, and the star formation histories (see Fig. 1121 . 

We simulate cluster and field galaxy samples containing 
1000 galaxies each. The results of the simulations in terms of 
observable quantities are summarized in Table |K| The mod- 
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els achieving relatively close values to the observed ones are 
marked with stars. The values of the scatter shown in the 
table represent the standard deviations for the simulated 
galaxy samples^, and the observational errors for the ob- 
served quantities. The standard deviations of the predicted 
quantities for the simulated field galaxies depend only on the 
assumed errors. As the table caption indicates, the standard 
deviation of AMs for the simulated field galaxies ai z = 0.4 
is 0.77, remarkably close to that of our observed field TFR 
(0.78; see Sec. 13. 11 . For the cluster galaxies, the scatter also 
depends on both the brightening of galaxies due to the burst 
of star formation and the fading of galaxies in the fade-out 
phase after the event. Nevertheless, assuming errors compa- 
rable to those of our observations, we find that the scatter 
in the observed properties of field and cluster galaxies are 
almost the same, and consistent with the observed ones. 
Since we are interested in relative differences between the 
cluster and the field galaxies, we ignore possible offsets be- 
tween the absolute colours of the models and the observed 
galaxies. For simplicity, we also ignore the fact that actual 
spiral galaxies show a wider spread in colour than our mod- 
els due to differences in morphology, inclination, extinction, 
stellar populations, metallicities, etc. We expect differential 
quantities to be more robust than absolute ones. 

As a first step, we explored the effects of varying some 
of the parameters on the predicted differences between the 
cluster and field model galaxies. We found that the depen- 
dence of the model results on the actual value of Zots is small 
(see Models 2 and 3). This is due to the fact that our ob- 
servations were designed to reach approximately the same 
absolute magnitude at each redshift. Thus using our fiducial 
value Zobs = 0.4 to model our complete dataset should be 
safe. The value of Atsp that best corresponds to our Subaru 
sample is not well defined by our galaxy selection procedure. 
Luckily, changing this parameter within very broad limits 
has almost negligible effect on the predicted differences be- 
tween field and cluster model galaxies (Models 4 and 5). 
Changing the magnitude limit of the observations has only 
a moderate effect (Model 6): going deeper in the luminosity 
function changes the results a little by including a few more 
galaxies in the fade-out phase in the galaxy samples. 

The main conclusion of exploring models 1-6 is that 
with the parameters explored there we do not get results 
which simultaneously agree with our observed changes in 
luminosity, colour and absorption-line galaxy fraction. A 
better match to our observations is achieved if we reduce 
the burst strength to 0.1 (Model 7). A truncated star 
formation history without a star-burst (Model 8) also pre- 
dicts close values to our observations in terms of magnitude 
and colour differences because cluster emission-line galax- 
ies spend a very short time in the emission phase, and thus 
do not have time to make their properties different from 
the field galaxies. However, the fraction of absorption-line 
galaxies ends up being too high for this model since the 
vast majority of the cluster galaxies would have stopped 
forming stars. Alternatively, our observed values are roughly 



The error in the mean values can thus be obtained by dividing 
the standard deviation by (N — 1)^/^, where N = 1000 is the 
number of galaxies in the simulations. 



achieved also if the burst timescale is as long as ~ 0.3 Gyr 
(Model 9). 

The magnitude and colour distributions from Models 1 
and 9 are over-plotted in Fig. |^ It seems as if Model 9 not 
only provides a better match to the average observed val- 
ues, but also to their distributions. However, our models 
tell us that there is a degeneracy between the burst mass 
fraction and the star-formation time scale. Models 7 and 9 
reproduce the observations almost equally well. Neverthe- 
less one interesting conclusion is that it is possible to have 
substantial starbursts in the cluster spirals without signifi- 
cantly changing the average colours and luminosities of the 
population of line-emitting spiral galaxies in a sample like 
ours. It is also clear that the scenario we propose would im- 
ply a significant increase in the fraction of spiral galaxies 
with absorption-line spectra when the field galaxies fall into 
the clusters, as observed. 

On the other hand, B05A found that the mean TFR 
residuals of the cluster galaxies are 0.7 mag brighter than 
that of the field galaxies, in agreement with the earlier 
findings of Milvang- Jensen et al. (2003). This is closer to 
the predictions of more massive or shorter timescale bursts 
(e.g.. Model 1) than our case. If we assume that the results 
from B05A can be represented by Model 1, the difference 
from the Subaru observations would be ~ 1.1a in terms of 
A^''^'"(AMb) (~ 1.6a when compared with B05A directly 
as opposed to comparing with the model). However, the dif- 
ference would be ~ 4.5cr in terms of A^.''^'"(_B — V), i.e., 
the model predicts significantly bluer cluster emission-line 
galaxies. Taken at face value, this implies a very low prob- 
ability that the discrepancies in the results from the VLT 
and Subaru data are due simply to a different sampling of 
the same parent population. Of course, this assumes that 
one can apply simple Gaussian statistics to this complex 
problem, which is far from clear. 

Thus, we cannot rule out that other factors not ac- 
counted for in our simulation could also play some role in 
the observed differences. The lower median redshift of the 
Subaru sample (0.39), as compared with the VLT one (0.52), 
would suggest that one would expect a stronger evolution- 
ary effect for the latter. However, B05A show that the effect 
is present for even their lower redshift clusters {z ~ 0.3). 
Cluster-to-cluster differences are, of course, another poten- 
tial candidate. 

The fact that we get a relatively weak statistical signif- 
icance for the inconsistency between the Subaru results and 
Model 1, or between the Subaru data and the VLT data of 
B05A in terms of A^.''^'"(AAfB) is due to the large uncertain- 
ties associated with the measured values. Indeed, according 
to our model simulations, one would need to determine reli- 
able rotation velocities for ~ 50 cluster spirals and a similar 
number of field ones to distinguish Model 1 (strong, short 
burst) from Model 9 (weak burst with slowly-declining star 
formation) at 3a level by using A^''^'"(AM_b) alone. Thus 
much larger samples than those presented here and in B05A 
are needed. 

5.2 Evolution of star formation rate with redshift 

The evolution of the B-band absolute magnitude with red- 
shift for the field spirals in our sample (Eq. Q and Fig. |SJ 
can be interpreted as due to the evolution of the star- 
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Table 5. Model predictions for different parameters: The values of the model parameters for the 
standard model arc discussed in the text. Only the parameters that change in each model are listed. The 
remaining parameters are as in the standard model. The models achieving values close to the observed ones 
are marked with stars. 


iVioucl ^ ana item 


Atejji 

(Gyr) 


(mag) 


A^'<=™(B - V)^ 
(mag) 


A='"'"(B - 1/)^ 
(mag) 


f K 

J c,abs 


(1) Standard model 

Zobs = 0.4, Atsp = 3 Gyr, 
Mii^ = M*(z), /t = 0.2, 
T = 0.1 Gyr 


n OK 

U.ZD 


— U.D4 ± u.oZ 


n 1 o ^ n 1 
— yj.LZ ± u. io 


U. io ± U. ly 




(2) ^ 0.2 

(3) 2^6, ^ 0.6^ 


0.26 
0.26 


-0.52 ±0.79 
-0.52 ±0.81 


-0.10 ±0.15 
-0.10 ±0.15 


0.26 ± 0.22 
0.15 ± 0.20 


53 
31 


(4) Msp 1 Gyr 

(5) Atsp 8 Gyr= 


0.26 
0.26 


-0.49 ±0.75 
-0.54 ±0.77 


-0.11 ±0.14 
-0.11 ±0.15 


0.14 ±0.17 
0.19 ± 0.21 


42 
46 


(6) Miin,^M*{2) + 2mag 


0.26 


-0.42 ±0.83 


-0.11 ±0.15 


0.33 ± 0.22 


76 


(7) A ^0.1 * 

(8) h O.O/ 


0.22 
0.002 


-0.12 ±0.77 
0.10 ±0.98 


-0.03 ±0.16 
0.02 ±0.19 


0.23 ± 0.19 
0.25 ± 0.18 


46 
99 


(9) T 0.30 Gyr * 


0.57 


-0.03 ±0.76 


0.00 ±0.16 


0.27 ±0.17 


30 


(10) Burst type : Exp.^ Const. 0.11 


-0.96 ±0.76 


-0.20 ±0.15 


0.13 ±0.20 


60 


Observed^ 




-0.18 ±0.33 


0.06 ±0.04 


0.20 ±0.05 


20-34 



Difference in TFR residuals between cluster and field galaxies, and the standard deviation for the cluster 
galaxies (see text). The average residual for the field galaxies and its standard deviation sA, z = 0.4 is 
AMgj = -0.50 ±0.77. 

^ Difi'erence in rest-frame B — V colour between the emission-line galaxies in the clusters and the field, and 
the standard deviation for the cluster galaxies. The average colour of the field galaxies and its standard 
deviation is (B -V)f = 0.34 ± 0.14. 

The same as b, but for the absorption-line galaxies. 

Fraction of absorption-line cluster galaxies. Note that the observed fraction shows the measured fraction 
of cluster absorption-line spirals minus the measured fraction of field absorption-line spirals. This is done 
because our models assume that all the field spirals have emission lines. 

^ The nominal ages of the model for the cluster galaxies can exceed the age of the universe, although such 
galaxies are in the fade-out phase and thus insensitive to the assumed age of the underlying population at 
the time of the burst. 

f No starburst (truncated star-formation). 

5 The observed mean values with 1 a errors except for /cats: where we give the accepted range. 



formation rate (SFR) of these galaxies. In this section, we 
parameterise the SFR evolution as a simple power law of 
the form SFR oc (1 + 2)", and estimate the value of a that 
best explains our result. In order to calculate the luminosity 
evolution in the B-band associated with the postulated SFR 
evolution we used Bruzual & Chariot (2003) models with a 
Chabrier (2003) initial mass function and solar metallicity. 
The results are shown in Fig. 1131 Changing the formation 
epoch of the galaxies in the range 1 < z < 5 hardly affects 
our estimated a for a < 2. Thus, we arbitrarily set the time 
at which spirals start forming stars at 20% of the Universe's 
present-day age {zf — 2.7 with the adopted 'concordance' 
cosmology). Note that if the SFR had remained constant 
with time {a — 0), spiral galaxies would get slightly fainter 
with redshift in the B band as the building up of their stellar 
populations with time slightly overcompensates the average 
ageing of their stellar populations. Since our observations 
indicate that spiral galaxies get brighter with redshift, their 
SFR must be higher at higher z, i.e., a > 0. 



To compare our observational results with the model 
predictions, and to increase the statistic accuracy, we com- 
pute the weighted-mean of the slopes in Eq. Q and in 
B05B (—1.0 ± 0.5), and interpret it as the representative 
of the evolution of field spiral galaxies at the median red- 
shift of the VLT and Subary field galaxy samples. We obtain 
AMb = —0.51 ± 0.22 mag at z = 0.48, which is shown in 
Fig. ll3l as the data-point with errorbars. Taken at face value, 
our analysis yields a = 1.7 ± 0.7. 

This value could be an overestimate of the rate of evolu- 
tion for two reasons. First, in our analysis we have assumed 
constant (solar) metallicity for the spiral galaxies. It is rea- 
sonable to expect that the average metallicity of field spirals 
might have been lower in the past, and, for a given stellar 
mass, a lower metallicity stellar population is expected to 
be brighter in B (see, e.g., Bruzual & Chariot 2003). Thus, 
part of the observed luminosity evolution may be caused by 
a change in the average metallicity of the stellar popula- 
tions of the galaxies. However, we expect this effect to be 
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Figure 13. The evolution of B-band absolute magnitude with 
redshift for galaxies with different star-formation histories. The 
lines show the model predictions with SFR oc (1 + z)" for differ- 
ent values of a and solar metallicity. The circle with the errorbar 
shows the slope derived from Eq. combined with the one de- 
rived by B05B placed at the median redshift of the combined field 
TFR sample. 

small because the measured evolution in the average metal- 
licity of field galaxies is rather modest (cf. Kobulnicky et al. 
2003; Lilly, CaroUo & Stockton 2003; Kobulnicky & Kew- 
ley 2004). At the median redshift of our field galaxy sample 
we expect the average metallicity (at a given galaxy mass) 
to be only ~ 0.07 dex lower than a.t z — (Kobulnicky 
& Kewley 2004), which would have negligible effect in our 
calculation. The second reason why we may have overesti- 
mated the rate of SFR evolution is that any unaccounted for 
selection effects probably mean that we preferentially select 
intrinsically brighter galaxies at high-z, thus increasing the 
apparent luminosity evolution. 

This rate of SFR evolution is substantially shallower 
than the rate of evolution of the SFR density of the universe 
indicated by studies based on the UV, Ha, far-infrared, and 
radio emission of galaxies, which all suggest a ~ 3-4 be- 
tween z = and 1 (e.g., Hopkins 2004). Hence, the rapid 
evolution in the cosmic SFR density is not driven by the 
evolution in the SFR of individual bright {Mb ^ Mg) spi- 
ral galaxies like the ones in our sample, in agreement with 
the conclusions of B05B. 



6 SUMMARY 

We have carried out MOS observation of 4 cluster fields us- 
ing the FOCAS spectrograph at Subaru, and obtained spec- 
tra of 103 cluster and field spiral galaxies with spectroscopic 
redshifts in the range 0.06 55 2 ^ 1.20. A total of 77 galaxies 
show emission lines. Of these, 33 galaxies yielded observed 
rotation curves of good enough quality to determine secure 
rotation velocities. Our sample reaches roughly Mg at each 
redshift. 

By comparing the rest-frame B-band TuUy-Fisher rela- 
tion (TFR) of our cluster and field emission-line spiral galax- 
ies, we found no measurable difference between the B-band 
luminosity at a given rotation velocity of both populations. 



This agrees with the conclusions of Ziegler et al. (2003), 
but disagrees with our own previous VLT results (Milvang- 
Jensen et al., 2003; Bamford et al. 2005A). We also find the 
rest-frame B — V colour of the cluster emission-line galaxies 
is marginally redder (by 0.06 ± 0.04) than that of the field 
galaxies, providing little indication that the cluster spirals 
are undergoing enhanced star formation. On the other hand, 
we find that the fraction of spiral galaxies with absorption- 
line spectra (i.e., no detectable emission lines) in the clusters 
is larger than that in the field by a factor ~ 3-5, in agree- 
ment with the results of Dressier et al. (1999). This implies 
that the cluster environment quenches star-formation on its 
spiral galaxies, leading, perhaps, to the formation of SOs. 

To evaluate the significance of our results, we carried 
out simulations of the effects that changes in the star- 
formation history of in-falling spiral galaxies would have 
on the observed properties of galaxy samples similar to 
ours. Following the encounter of the in-falling galaxies with 
the cluster environment, our models explore different star- 
formation scenarios, including truncation of the galaxies' 
star-formation on different timescales, possibly preceded by 
bursts of different strengths. It turns out that even quite 
drastic changes in the SFR of the galaxies may have quite 
modest effect on the average luminosity and colour of a sam- 
ple of bright cluster spiral galaxies selected to have ongoing 
star-formation (i.e., having the strong emission lines needed 
for rotation-curve measurement). However, the fraction of 
spiral galaxies with absorption-line spectra (i.e., without 
current star-formation) is very sensitive to the effect of star- 
formation truncation. Our Subaru observations favour mod- 
els with relatively mild or absent initial star-bursts, and rel- 
atively long star-formation timescales, while the VLT results 
of Bamford et al. favour a more massive initial star-burst and 
a shorter time-scale. However, we estimate the probability 
that the observed difference between our Subaru results and 
the VLT ones of Bamford et al. arises from "unlucky" statis- 
tic sampling of galaxy populations with intrinsically similar 
properties. We find that both TuUy-Fisher relation results 
are only different at the ~ la level. This is due to the rela- 
tively small sample sizes and the large uncertainties in the 
determination of the individual TFR offsets. To definitively 
rule out the presence of a starburst before the star-formation 
cessation would require samples of ~ 50-100 field and clus- 
ter spirals with reasonable determinations of their rotation 
velocity. 

Finally, we find that the rest-frame absolute B-band 
magnitude (at a fixed rotation velocity) of the field galax- 
ies in our sample shows an evolution of —1.30 ± 1.04 mag 
per unit redshift. By statistically-combining our luminos- 
ity evolution estimate with that of Bamford et al. (2005B; 
— 1.0 ±0.5 mag per unit redshift), and interpreting it as due 
to the increase with redshift of the SFR of the galaxies, we 
estimate that SFR cx {1 + z)^'"''^'^''^ for our field spirals. This 
indicates that the average SFR of bright {Mb < Mg) disk 
galaxies evolves more slowly than the universal SFR density, 
suggesting that the evolution of the global SFR evolution is 
not dominated by bright star-forming disk galaxies, in agree- 
ment with previous studies. 
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